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ABSTRACT 

It  is  shown  that  for  an  arbitrary  strictly  increasing  knot  sequence 
00 

t = (t  ) and  for  every  i,  there  exists  exactly  one  fundamental  spline 
— ' \ — 00 

L (i.e.,  L (t  ) = 6 , all  i),  of  order  2r  whose  r-th  derivative  is 
i ’ i i if 

square  integrable.  Further,  l|F  (x)  is  shown  to  decay  exponentially  as  x 
moves  away  from  t.,  at  a rate  which  can  be  bounded  in  terms  of  r alone 
This  allows  one  to  bound  odd-degree  spline  interpolation  at  knots  on 
bounded  functions  in  terms  of  the  global  mesh  ratio  M :=  sup  At  /At  . 

i 1 » J l 

A very  nice  result  of  Demko's  concerning  the  exponential  decay  away 
from  the  diagonal  of  the  inverse  of  a band  matrix  is  slightly  refined  and 
generalized  to  (bi)infinite  matrices. 
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Odd-degree  spline  Interpolation  at  a bllnflnlte  knot  sequence 

Carl  de  Boor* 

1.  Introduction.  Let  i :=  (t^)^  be  * bllnflnlte,  strictly  Increa- 
sing sequence,  set 


%.  = = tl 

1-*!® 


let  k = 2r  be  a positive,  even  Integer,  and  denote  by  + the  collect- 

* 

0 “ 

Ion  of  spline  functions  of  order  k (or,  of  degree  < k)  with  knot  se- 
quence £.  Explicitly,  ^ consists  of  exactly  those  k-2  times  contin- 
uously differentiable  functions  on 


I **  V 

which,  on  each  Interval  (t^ , coincide  with  some  polynomial  of 
degree  < k,  l.e. , 


*i,s  ■*  'k.j'10'1'2  ~ 1 ■ (t-»-  V • 

We  are  particularly  Interested  In  bounded  splines 

“Vi  !=  ViAB(I)* 

l.e.,  In  splines  s for  which 


IlsIL  :=  »up  | a ( t)  | 

® t«I 

Is  finite.  It  Is  obvious  that  the  restriction  map 
Ri  S ♦k.i-^5  != 

oarrles  m$^  ^ Into  the  space  m(Z)  of  bounded,  bllnflnlte  sequences. 
Ws  are  Interested  In  inverting  this  map,  l.e..  In  Interpolation.  We 
consider  the 

Bounded  Interpolation  Problem:  To  construct,  for  given  a tm(Z)  . 
some  s c ml^  ^ for  which  s | ^ = a . 

* Sponsored  by  the  United  States  Army  under  Contract  DAAG29-75-C-0024 
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we  will  say  that  the  3.I.P.  lo  correct  (for  the  given  knot  sequence 

|)  If  It  hasexactly  one  solution  for  every  oen(Z). 

We  consider  under  what  conditions  on  1 the  B.I.P.  Is  correct.  We 

also  discuss  the  continuity  properties  of  the  map  at— »s  In  case  the 

a 

B.I.P.  Is  correct.  We  establish  the  following  theorem. 

Thejjrem_l_.  If  the  global  mesh  ratio 
M.  :=  sup  At. /At . . 

5 i,J  1 3 ' . 

Is  finite,  then  I = (•<»,  a>) , artf.  R^  maps  mf^-  faithfully  onto  m(Z)  , 
l.e..  for  every  bounded,  bllnflnlte  sequence  a,  there  eilsts  one  and  only 
3M  ...bg.U&itefl  SPllae  for  which  s(x(t1)  = o1#  all  i.  Moreover. 

(i.l)  IsX  < const  Hall^,  all  afem(Z)  , 

with  const  depending  only  on  k and  M.  . 

1 

We  note  In  passing  the  following  Immediate  corollary. 

o 

Corollary , Denote  by  0[a,b]  the  space  of  continuous  (b-a)-pcrlodlc 
functions  on  R.  Given  I :=  (T^jj  with  a = TQ  < ...  < rQ  = b,  let  i = 
(t^)^  be  Its  "(b-a)-perlodlc  extension".  l.e.. 

j • — J (b— a)  for  1=1 ....  .n  and  all  J a2Z« 

B*no*e  by  tilt.  (b-a)-perlodlc  functions  In  1^.  Then  (as  Is  well 

known)  , £or  every  f tc[a,b]  , there  exists  exactly  one  sf  e which 

With  f JLt  r0 , tj_ , ...,  rn.  further,  for  some  const  depending  only 
ga_.the  global  mesh  ratio  Mj.  = max1>3  A^/ATj, 

l»fl®  - const  If f , ftl*.  f «C[a,b]. 

Indeed,  if  agrees  with  feO[a,b]  at  $,  then  so  does  its 

translate  s^(*  - (b-a))  which  is  also  in  and  thersfors  must  equal 

•e,  by  the  uniqueness  of  the  interpolating  spline.  This  shows  that  sf  la 
the  interpolating  spline  in  for  f,  and  so  Rs,!  < const  ||f||  from  (1.1). 

for  the  case  of  uniform  i,  i = Z say,  the  problem  of  bounded  In- 


terpolatlon  has  been  solved  some  time  ago  by  <Ju.  Subbotln  [17].  In  this 
case,  the  Interpolation  conditions  *a|j  = a ®stabllsh  a one-to-one 
and  continuous  correspondence  between  bounded  splines  and  bounded  se- 
quences. Subbotln  came  upon  the  Interpolating  spline  as  a solution  of 
the  extremum  problem  of  finding  a function  s with  s|^  = o and  smallest 
possible  (k-l)st  derivative,  measured  In  the  supremum  norm.  Later, 

I.J.  Schoenberg  investigated  the  B.I.P.  once  more,  this  time  as  a spe- 
cial case  of  cardinal  spline  interpolation. to  sequences  a which  do  not 
grow  too  fast  at  Infinity  [15],  [16], 

Little  is  known  for  more  general  knot  sequences.  The  simplest  case, 
It  * 2,  of  piecewise  linear  interpolation  Is,  of  course,  trivial.  The 
next  simplest  case,  k = 4,  of  cubic  spline  Interpolation  has  been  In- 
vestigated in  [6]  where  the  above  theorem  can  be  found  for  this  case. 

The  basic  tool  of  the  investigation  In  [6]  is  the  exponential  de- 
cay or  growth  of  nullspllnes.  Nullspllnes  are  therefore  the  topic  of 
Section  2 of  this  paper,  if  only  to  admit  defeat  in  the  attempt  to  ge- 
neralize the  approach  of  [6].  We  are  more  successful,  In  Section  3,  In 
identifying,  for  each  knot  sequence  £ and  each  1,  a particular  funda- 
mental spline  L^ , i.e.,  a spline  with  L^Ct^)  = G^,  which  must  figure 
In  the  solution  of  the  B.I.P. , if  there  is  one  at  all  (see  Lemmas  1 and 
2).  The  argument  Is  based  on  an  idea  of  Douglas,  Dupont  and  Wahlbin  [12] 
as  used  In  [7]  and  further  clarified,  simplified  and  extended  by  S.Demko 
[lO].  It  Is  also  shown  (in  Lemma  3 and  its  corollary)  that  the  r-th  de- 
rivative of  a nontrivial  nullspllne  must  Increase  exponentially  In  at 
least  one  direction.  The  exponential  decay  of  the  fundamental  spline 
is  used  in  Section  4 to  prove  Theorem  1.  That  section  also  contains 
a proof  of  the  fact  (Theorem  4)  that  the  B.I.P.  is  solvable  in  terms  of 
exponentially  decaying  fundamental  splines,  if  it  is  oorrect  at  all. 

This  fact  is  closely  oonneotad  with  S.Demko's  results  [lO]. 
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2.  Nullspllnes  and  fundamental  splines.  It  la  clear  that  the  prob- 
lem of  finding,  for  an  arbitrary  given  bllnflnlte  aequence  a,  eome 
spline  t for  which  sj^  = a,  always  has  solutions.  In  other  words. 

It  Is  clear  that  R^  maps  ^ onto  I?.  To  see  this,  start  with  a poly- 
nomial pQ  of  order  k which  satisfies  P0(tQ)  = aQ,  p^t^)  = , and  set 

s a pQ  on  [tQ,  t^].  Now  suppose  that  we  have  s already  determined  on 
some  Interval  [t^.t^]  and  let  Pj_^  be  the  polynomial  which  coincides 
with  s on  [tj_^ftj].  Then 

Pj(t)  :=  P^U)  ♦ 

is  the  unique  polynomial  of  order  k which  takes  on  the  value  Bj+j.  at 
tj+1  and  agrees  with  p^_1  (k-l)-fold  at  The  definition 

■ = Pj  °n  Ctj»tj+l^ 

therefore  provides  an  extension  of  s to  C "t^  * and.  In  fact,  the 

only  one  possible.  The  extension  to  is  found  analogously. 

In  this  way,  we  find  a solution  inductively. 

The  argument  shows  that  we  can  freely  choose  the  interpolating 
spline  on  the  interval  [t^,t^]  from  the  k-2  dimensional  linear  mani- 
fold 

\p«JPk  : P(t0)=a0*  P(V=al^ 

and  that,  with  this  choice,  the  interpolating  spline  is  otherwise  uni- 
quely determined.  In  particular,  the  set  of  solutions  for  o = 0,  l.«., 

the  kernel  or  nullspace  of  the  restriction  map  R+ , is  a k-2  dimensio- 

i 

nal  linear  space,  whose  elements  we  call  nullspllnea.  In  other  words, 
nullspllnes  are  splines  which  vanish  at  all  their  knots. 

The  difficulty  with  the  B.I.F.  is  therefore  not  the  construction 
of  some  interpolating  spline.  Rather,  the  problem  is  interesting  be- 
cause we  require  an  interpolating  spline  with  certain  additional  charw 
acteristics  or  "side  conditions" , viz.  that  it  be  bounded.  Nullspllnes 


4 


can  be  made  to  play  a major  role  In  the  analysis  of  this  problem. 

?or  Instance,  the  question  of  how  many  bounded  solutions  there 
are  Is  equivalent  to  the  question  of  how  many  bounded  nullspllnes  there 
are.  More  Interestingly,  a well  known  approach  to  the  construction  of 
lnterpolants  consists  In  trying  to  solve  first  the  special  problem  of 
finding,  for  each  1,  a fundamental  spline,  i.e.,  a spline  e ^ for 
which 


L^(tj)  = ®l«  J * ®X1  3*  , 

Such  a spline  consists  (more  or  less)  of  two  nullspllnes  Joined  together 
smoothly  at  t^.  Therefore,  If  one  could  prove  that  both  nullspllnes  de- 
cay exponentially  away  from  t^ , l.e., 


I|Ll”(tj'tj4l)  ~ C°n3tk  Xll"Jl 
at  a rate  X i [0,1)  whloh  Is  Independent  of  1, 

the  series 


, all  J , 

then  It  would  follow  that 


® 

(2.1)  a :=  E a.L. 

° 1=-®  1 1 

converges  uniformly  on  compact  subsets  of  I and  gives  a solution  sQ  to 
the  G.I.F..  In  fact,  sQ  theadepends  continuously  on  a,  l.e., 

- constk,X  » #11  ®tm(2) 

for  some  const^  ^ which  does  not  depend  on  a. 

The  hope  for  such  exponentially  decaying  fundamental  functions  is 
really  not  that  farfetched.  Suoh  functions  form  the  basis  for  Schoen- 
berg's analysis  In  the  case  of  equidistant  knots,  and  they  occur  Impli- 
citly already  In  Subbotln's  work.  Further,  a very  nice  result  of  S.Demko 
[10]  .3  be  elaborated  upon  In  the  next  section  (see  also  C.Chul's  talk 
at  this  conference)  shows  that  the  bounded  spline  lnterpolant  sQ  to 
bounded  data  a Is  neoessarlly  of  the  form  (2.1)  with  exponentially  de- 
caying In  case  sa  depends  continuously  on  a. 

In  a rather  similar  way,  nullspllnss  also  occur  In  the  discussion 


-S- 


of  Interpolation  error.  If  f 13  sufficiently  smooth,  end  sf  Is  Its 
spline  lnterpolant,  i.e.,  sf|t  = f|t  * then  one  gets,  formally  at 


flret,  that 


(2.2) 


- sf(t)  = / °°  K(t,s)  f( 


(s)  ds  . 


Here,  the  Peano  kernel  K(t,*)  Is  a spline  function  of  order  k with 
knots  t and  an  additional  knot  at  the  point  t,  and  vanishes  at  all  the 
knots  i.  Hence,  K(t,0  Is  again  a function  put  together  from  two  null- 
splines.  The  exponential  decay  of  these  two  nullspllnes  away  from  t Is 
desirable  here,  since  only  with  such  a decay  can  (2.2)  actually  be  veri- 
fied for  Interesting  functions  f.  But,  I won't  say  anything  more  about 
this  here. 

Based  on  my  experience  with  [6],  I had  at  one  time  considerable 
hope  that  the  exponential  decay  of  nullspllnes  could  be  proved  with  the 
help  of  the  following  considerations.  A nullspline  s c,  1 + Is  determined 
on  the  interval  Ct^,t1+^]  as  soon  as  one  knows  the  vector 

ai  :=  («,(t1) s(k-2)(t1)/(k-2)*) 

since  one  knows  that  s(t^)  = s(t^+1)  = 0.  One  can  therefore  compute  s1+1 
from  Sj^  In  a linear  manner.  Specifically, 

®1«-1  = a*  * 

with  A(h)  the  matrix  of  the  form 

A(h)  :=  diagd.h"1,...  ,h"k+3)  A diag(l  ,h,. . . ,hk_3) 

and  A = A ( 1 ) the  matrix 

* ■■  UkiKJ)n:U- 

This  means  that  A(h)  has  many  nice  properties.  For  instance,  A"*(h)  = 
A(-h) , and  A(h)  is  an  oscillation  matrix  in  the  9ense  of  Gantmacher  and 


Krcln. 


In  the  special  cubic  case,  k = 4,  A(h)  has  the  simple  form 


-6- 


A(h)  = 


2 


h 


3/h  2 

A 

and  allows  therefore  the  conclusion  that  s^  grows  exponentially  either 
for  Increasing  or  else  for  decreasing  Index  1 , at  a rate  of  at  least  2. 
This  observation  goes  back  to  a paper  by  Blrkhoff  and  the  author  [l]. 

The  transformation  1(h)  has  been  studied  In  much  detail  in  the 
oase  of  equidistant  knots  In  a paper  by  Schoenberg  and  the  author  [8], 
and  also.  In  more  generality,  by  C.  Micchelll  [14],  But,  such  exponent- 
ial decay  or  growth  for  nullspllnes  on  an  arbitrary  knot  sequence  has 
so  far  not  been  proved.  S.  Frledland  and  C.  Micchelll  [13]  have  obtained 
from  such  considerations  results  concerning  the  maximal  allowable  local 
mesh  ratio 

m*  :=  sup  At, /At . . 

^ 11-31=1  1 3 

3.  Exponential  decay  of  the  r-th  derivative  of  fundamental  splines 
and  nullspllnes  of  order  k = 2r.  We  base  the  arguments  In  this  section 
on  the  best  approximation  property  of  spline  Interpolation.  To  recall, 
the  r-th  divided  difference  of  a sufficiently  differentiable  function 
f at  the  points  t^,  . ..,  t^+r  can  be  represented  by 

[t1#...,t1+r]f  = jM1(t)  f(r)(t)  dt/rJ 

with  M.  = Ms  a B-spline  of  order  r, 

* ififi 

M1(t)  :=  r[tlf...,t1+r](.  - tjJ"1, 

normalized  to  have  unit  Integral.  Further,  (s^r*  : s c $2r  = $r  t 

while,  by  a theorem  of  Curry  and  Schoenberg  [9], 

*r,i  = { VlMl  * on  1 • 

where  we  take  the  bllnflnlte  sum  polntwl se , l.e., 

(EiP1M1)(t)  s=  Ej^M^t),  all  tfeK. 

This  makes  good  sense  since 


T 


M^(t)  > 0 with  strict  Inequality  Iff  t^  < t < t^+Jl  . 

Lemna  1.  Let  :=-(L€$2r  t : L(tj)  = *11  j\«  Then  £,^  has 

exactly  one  element  In  common  with  g.jr^  (I) . We  denote  this  element  by 

L1 

and  call  It  the  1-th  fundamental  spline  for  the  knot  sequence  J.  Further. 
with  the  abbreviations 

(3.1)  E :=  sup  At.,  h :=  Inf  At.  > 

3 J 3 J 

we  have 

(3.2)  llL{r)|l2  < constr  ^2/hT 
for  some  constant  constr  depending  only  on  r. 

Proof^.  We  first  prove  that  contains  at  most  one  element  In 
X*r)(I)  = (fftCr_1(I)  : f(r”1)  abs.cont. , f (r)  c X2(I)\ . Since 
= ker  Hj,  it  Is  sufficient  to  prove  that  the  only  null  spline  In  I»2r^  18 
the  trivial  nullspline.  For  this,  let  sCker  R^.  (I) » Then,  by  the 

Introductory  remarks  for  this  section, 

s^  = RjPjMj  for  some  pci?,  s^cL.,,  and  f M^s^  = 0 for  all  J. 

But,  by  a theorem  In  [3],  there  exists  a positive  constant  Dy  which  de- 
pends only  on  r so  that,  for  1 < p < ®,  and  for  all  rel? • 

(3.3)  d;1  Hr llp  < II E-jY J ( ( 1 3 +r- 1 3 ) /r ) 1"1/pH 3II p < llrllp  . 

Here,  ||rllp  :=  ( Z ^ | r I P ) 1/^P  > while,  for  f on  I,  ||fHp  : = (/x  |f  lp)1//p.  This 
shows  that  the  sequence  (Mj)  given  by 

(3.4)  Mj  :=  ((tj+r-*j>/r)1/a  «j.  all  3. 

A 

Is  a Schauder  basis  for  $p  jr\X2*  Therefore,  EjYjMj  converges  I«2  to 
the  spline  function  In  &2  It  represents.  But  this  means  that  our  part- 
icular spline  s(r)  is  In  the  3L2-span  of  (M^) , yet  orthogonal  to  every 
one  of  the  , which  means  that  e^r^  vanishes  identically.  But  then, 
since  s vanishes  more  than  r times, s Itself  must  vanish  Identically. 

-8- 
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Next,  we  prove  that  contains  at  least  one  element  In  l^r^(I). 
For  this,  we  recall  from  [5]  that  there  exists,  for  any  given  a«I?, 
a function  g which  Is  locally  In  JL^and  satisfies  gjt  = a,  and  whose 
r-th  derivative  satisfies 


J«)2)l/2. 


(3.5)  ilg  r Il2  < Dr(l  (^j+r-tj)(Ctj,...,tJ+r]a)2)1/2, 

with  Dr  the  same  constant  mentioned  In  (3.3).  Here,  the  number 
[ t j ,. . . ,t j+r]a  stands  for  the  rth  divided  difference  at  the  points  t^, 
...»  'tj+r  of  any  function  f for  which  f|^  = a.  In  this  way,  we  obtain 
for  the  specific  sequence  a = (6^_  j)  a function  gfeX^  for  which 

gftj.)  = 6^  j,  all  J, 

while  llg^r^ll2  is  bounded  by  the  right  side  of  (3.5).  Note  that,  for 
the  specific  sequence  a = (6^  ^),  this  bound  becomes 


- Br(JJ.JVr-*J)Cl/^}(Vtn>l2)1/2 


< constr  (h)^ /hr  . 

Now  let  g be  any  element  In  so  that  g^  is  the  I^-approxlmatlon 

to  g^  from  . f\TL0.  This  makes  sense  since  (3.3)  Insures  that 
$p  ^ f\ ~2  “ closed  5'uuSpaC c oT  X 2 ( I ) • Then 


t.+I>]g  = = /Mjg(rVr!  = [tj,...,t 

J,  while  Bg  r^||2  < || g ^ || o«  But  this  means  that,  for  an  appri 


••*wJ+rJ6* 

appropriate 


polynomial  p of  order  r. 


(s  ♦ pHt^)  = = 6i-j»  a11 

while  still  ||  (g  ♦ p ) ^ r ^ II 2 < ||g^ll2  < const^1/2/^1.  This  shows  that 

L :=  6 + P is  a function  of  the  desired  kind.  ||| 

We  continue  to  use  the  inequality  (3.3)  and  the  abbreviation 
Mj  = ( (tj+r-t  j)/r)1//2Mj , and  come  now  to  what  I consider  to  be  the 
main  point  of  this  paper. 


( r) 

Lemma  2.  If  0 Is  the  sequence  of  coefficients  for  L£  ' with  respect 
to  the  basis  (Mj  for  $_  _ , l.e..  If  L,r^  = L..8.M...  and 

j r ,j  i j j j 


,(n) 


0 , |j-l|  < n 

LP..  1 3-1 1 > n 


• n— 0 ,1 ,2 , . . . , 


then  there  exist  constr  and  Xp  €.  [ 0 , 1 ) depending  only  on  r 30  that 


(3.6) 


IIP 


(n), 


2 < constr  ||p||2  A°  . n=0,l,2 


The  Inequalities  (3.3)  allow  us  to  conclude  from  Lemma  2 the  ex- 
ponentlal  decay  of  L^  in  the  following  form. 

Corollary.  For  some  constr>  and  some  >y£[0,l)  depending  only  on  r, 
and  for  all  1 and  n, 

(r>. 


II Ll. 1 ^ II 2 , ( t_CD  » "ti_ n 1 + ^l^StUi+n-V  “ constr  llI'jLr;ll2  Ay 
Proof  of  Lemma  2.  Let 


(r). 


A :=  (fMjMj) 


be  the  Gram  matrix  for  our  appropriately  normalized  B-splins  basis  of 
$r  A proof  of  the  lemma  can  be  obtained  directly  from  the  fact  that 
the  elements  of  the  Inverse  matrix  for  A decay  exponentially  away  from 
the  diagonal  at  a rate  which  can  be  bounded  In  terms  of  r and  Indepen- 
dently of  £.  This  Is  proved  In  [7]  with  the  aid  of  a nice  Inequality 
due  to  Douglas,  Hxpont  and  Wahlbln  [12],  But,  between  the  time  I proved 
Lemma  2 this  way  and  the  delivery  of  this  talk,  S.  Demko  wrote  a paper 
[10]  In  which  he  demonstrated  that  such  arguments  use  actually  very 
little  specific  Information  about  splines.  Using  the  inequality  of  Doug- 
las, Dupont  and  Vfahlbln,  he  proved  the  following  nice 

Theorem  (3.  Demko).  Let  A :=  (a^)  be  an  Invertible  band  matrix 
(of  finite  order).  Explicitly,  assume  that,  for  some  m.  a^  = 0 when- 
ever 1 1—  J | > m.  and  that,  for  some  positive  K and  ?,  and  some  pe[l^»], 

5 Nip  < ii Axil p < riiznpt 
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all  x 


Then  the  entries  of  the  Inverse  A-1  =:  (b^j)  satisfy 
I bi  j | < const  X'1^'  • all  l.J. 

for  some  const  and  some  X *■  [ 0 ,1 ) which  depend  only  on  mt  p,  K and  K.  In 
particular,  these  constants  do  not  depend  on  the  order  of  the  matrix  A. 

The  Interested  reader  will  have  no  difficulty  In  proving  this  theo- 
rem after  a study  of  the  following  proof  of  Lemma  2,  a proof  which  makes 
essential  use  of  Demko's  ideas,  even  though  the  Inequality  of  Douglas, 
Dupont  and  Wahlbln  fails  to  make  an  explicit  appearance.  In  the  bargain, 
the  reader  will  thereby  obtain  explicit  estimates  for  const  and  X (which 
Demko  did  not  bother  to  compute). 

We  note  that  the  specifio  matrix  A = (Jm^M^)  16  a band  matrix,  of 
band  width  m = r-1  In  the  sense  that  = 0 for  |i-j|  > r-1.  Also, 

we  conclude  from  (3.3)  that  the  sequence- to- sequence  transformation 

a t— *Aa 

induces  a linear  map  on  -£2(Z)  to  /2(Z)  which  we  also  call  A and  wbloh 
1b  bounded  and  boundedly  lnv«£tlble.  Specifically,  one  obtains  from 
(3.3)  that 

0.7)  k:=  nA||2||A-1n2  < Dr2  . 

Here,  ||b||2  :=  sup  {||BaO g/lkllg  s a 6 ^(2Z)\ , as  usual. 

We  now  claim  that, 

(3.8)  for  all  n > 2r,  ||p(n)0^  < (K2/(l+K?))  l|P(n"2n)ll| 

which,  with  the  ^-independent  estimate  (3.7)  for  K,  establishes  the 
lemma  (with  X*  <U/(l+K2)l/2)l/2m  ). 

7or  the  proof  of  (3.8),  we  oonslder  without  loss  of  generality  only 
the  specific  function  Lq.  We  note  that 

(A0)^  - /M1Lj)r)  = r!( (t1+r-t^)/r)  [t^ ,.. • »t^+r]Lg 

= 0 unless  < tQ  < t1+r  . 
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Therefore , 

(3.9)  supp  Ap  c [-r,0]  , 

where,  for  any  bllnflnlte  sequence  a,  we  use  the  abbreviation 


supp  a :=  (i«2  s 
We  claim  that,  for  n > m, 

(3.10)  supp  Ap^  c (-n-m,  n+m)  \ (-n+m,  n-m)  . 

Indeed,  supp  (P^n*  - P)  c.  (-n,n) , hence  supp  A(p^n*  - P)  £ (-n-m,  n+m) 
whfch  also  contains  supp  Ap  = [-r,0],  therefore 

supp  Ap^n^  c (-n-m,  n+m)  . 

On  the  other  hand,  supp  P^  cZM-n»  n)  , therefore  also 
supp  Ap^  c 2N.(-n+m,  n-m}  . 

It  follows  from  (3.10)  that,  for  n > 2r, 

(3.11)  eupp  Ap^a*  0 supp  Ap^0"2®^  = / , 

therefore 

||AP(n>fl|  < ||AP(n)|||  ♦ HAP  (&-2B)H|  = l|A(P(n)-P(n‘2m>)l! 

But  then 

|U-1Hj1l|8(n,ll2  « 01f(nllla  < Il<(<nl-B(,,'5"l>ll2 

< Hl2»8(n)  - 8<n-an,lla  - 

l.e. , 

Dp(n)B|  < k2Hp(q)  - p(n-2a)ll| 

= R2(|p(n"2m)||  - flp(n)l|) 

which  proves  our  earlier  olaim  (3.8).  Ill 

It  Is  dear  that  the  argument  provides  the  exponential  decay  of 


<v  OJ 


the  form  (3.6)  and  with  X<  (k/(1+  k2)1^2)1^213  for  any  sequence  p In 

Jl^i Z)  for  which  Ap  has  finite  support.  In  particular,  one  obtains  such 

exponential  decay  for  the  sequence  for  which  Ay^  = (6^_j),  l.e.» 

for  the  1-th  row  of  the  matrix  inverse  of  A.  Further,  It  Is  clear  that 

(3.11)  Implies  ||Ap(n)||P  < ||Ap(n)||J  ♦ ||Ap(n_2m) ||^  for  any  1 < p <», 

hence,  the  argument  carries  at  once  from  ^(Z)  over  to  any  ^(Z)  with 

1 < p <«,  Denko  obtains  such  exponential  decay  also  for  p = ® by  con- 

T 

sidering  the  transposed  matrix  A for  which  then  automatically 
BAT||1B(AT)-1||1  = IIAII^IIA-1^ 

due  to  the  finite  order  of  the  matrix  he  considers.  This  switch  requires 
a word  or  two  in  the  Infinite  case*  as  follows.  As  one  easily  checks,  If 
a (bl) Infinite  matrix  (a^)  gives  rise  to  a bounded  linear  map  A on 
then  its  transpose  gives  a bounded  linear  map  B on  and  the  adjoint 
of  B is  then  neoessarlly  A Itself.  This  Implies  that,  If  a matrix  (a^) 
gives  rise  to  a bounded  linear  map  on  ^ which  Is  boundedly  invertible, 
then  Its  Inverse  can  also  be  represented  by  a matrix,  viz.  the  transpose 
of  the  matrix  which  represents  the  Inverse  of  the  linear  map  on  given 
by  the  transpose  of  (a^).  Of  course,  exponential  decay  away  from  the 
diagonal  is  unchanged  when  going  over  to  the  transpose. 

These  comments  establish  the  following 

Theorem  2.  Let  M be  a finite.  Infinite  or  bllnflnlte  "Interval"  In 
Z,  let  1 < p < ® , and  let  q :=  min  $p,  p/(p-l)\  . Let  (mi^  1 , JmM  be  a 
matrix  with  band  width  m :=  sup  { |1- J | s a^  r and  assume  that  (*j  j) 
Induces  a bounded  linear  map  A on  £p(M).  If  A Is  boundedly  Invertible, 
then  A"1  Is  also  given  by  a matrix,  (b^j)  say,  and 

l*ljl  * const  , all  1.1. 


I *• 


with 


«4 


X :=  (K/(l*it<l)1/<l)l/2m.  const  < !U-1||p/X2m.  K:=  ||A!!pBA“1Hp  . 

We  add  one  more  remark.  With  the  appropriate  Interpretation  of 
"handedness" , the  above  argument  carries  through  even  for  matrices  which 
are  not  banded  In  the  straightforward  sense.  As  a typical  example,  con- 
sider the  Gram  matrix  for  a local  support  basis  of  some  space  of  func- 
tions of  several  variables.  Then,  there  Is  no  ordering  of  that  basis  for 
which  the  corresponding  Gram  matrix  Is  appropriately  banded.  Eut,  If  we 
follow  the  geometry  of  the  underlying  problem  and  think  of  the  Gran 
"matrix"  as  acting  on  functions  on  some  multidimensional  Index  3et  M 
having  an  appropriate  metric  |.|  (Instead  of  on  TL)  , then  the  statement 
and  the  proof  of  Theorem  2 go  through  otherwise  unchanged.  We  do  not 
pursue  this  point  here  further,  but  alert  the  reader  to  Descloux's  fine 
paper  [ll]  In  which  such  considerations  can  be  uncovered  once  one  knows 
what  to  look  for. 

We  finish  this  section  with  the  observation  that  the  r-th  deriva- 
tive of  a nontrivial  nullspllne  must  Increase  exponentially  In  at  least 
one  direction.  The  argument  Is  rather  similar  to  the  proof  of  Lemma  2. 

We  continue  to  denote  by  A the  specific  matrix  (Jm^Mj)  and  recall 


(3.7)  K :=  llA|l2BA-1||2  < Dr2  . 

Lemma  3.  If  Is  the  r-th  derivative  of  a nullspllne  In  $2r,i 

and  1 < j are  arbitrary  Indices,  then 

(1  ♦ K2)  Z |0j2  < K2  ^ |pj2  . 

v=l  w=l-2m 


Proof.  Define  p',  0"  by 

f P*  » 

Pi  s=  ^ » 

(0  , otherwise 
•o  that  the  inequality  to  be  proved  reads 


0"  : 


- f *-• 1 

l 0 , o 


-2m<v<J*2m 


otherwise 


(1  ♦ K2)flp'|f  < K2  |0"lf 
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(3.12) 


Ve  have 


while 


■upp  Ap'  c [i-m,  Jt-m] 


■upp  (P  - p")  c Z , J+2o]  , 

therefore,  with  Ap  = 0, 

■upp  AP"  = aupp  A (p  - p")  C2\[U,  j+m]  c Z Ssupp  Ap ' . 
Consequently, 

flA-1!! flip'll 2 < ||AP’|I2  < tlA(p'  - p")||2  < ||a||2!1p'  - p"i2  , 
or,  with  K=  ||A||2||A”1||2, 

np'll|  < K2  IIP'  - P"||f  = K2  (l|P"l|f  - llp'llf) 
which  implies  (3.12).  | | | 


Corollary.  Let  I^PjM^  be  the  r-th  derivative  of  a null  spline  s in 
♦2r,i  >n<i  9et 

al  5=  E , IP.I2,  all  jtZ  , 

3 2mJ<i<2m( J+l)  1 

*ith  m : = iwi , as  before.  Then 

^•13)  Z a^  < *2  (a.  ♦ a.),  for  all  1 < j. 

i<v<J  * J “ 


Therefore,  for  all  ft,  and  either  for  all  i > 


1 < ft.. 


a^  > const^. 


with 


eonst^  :=  y a^/(A2A> 
and 

At=  (1  ♦ K2)/*2  > 1 . 


Proof.  Assertion  (3.13)  follows  at  onos  from  the  lemma.  The  second 
osssrtlon  of  the  oorollary  is  loss  obvious,  for  its  proof,  assume  with- 
out loss  that  ft-  s 0,  Prom  (3.13)* 


-I  s- 


Therefore , 
(3.14) 


A ao  i C aw , 1=1  *2 ,3 > •••  • 
-l<v<l 


t\  o < £ ao • 1~1  *2,3,  ... 

-1<»<1 


o :=  a^/A  . 

X 2 

Let  now  constQ  = ^ c/k  , as  defined  above,  and  assume  that  the  lnequal- 


a^  > constpA^I 
is  violated  for  some  1 > O while  also 


(3.15) 


a_j  < cons') 


for  some  positive  J which  we  assume  without  loss  of  generality  to  be  no 
less  than  1.  Then,  we  oan  also  assume  that  J Is  the  smallest  Index  > 1 
for  which  (3.15)  holds.  We  obtain  from  (3.13)  that 

(3.16)  £ »v  < fc.2(a  , ♦ a.)  < K2constJ  ♦ A1)  = \ c(A^  ♦ A1). 

-i<v<l  ~3  0 

On  the  other  hand,  by  (3.13)  and  by  the  choice  of  J, 

£ au  ™ £ **  ♦ £ a>, 

-3<v<i  -3<VS-1  -1<V<1 

» oonstjl  A^-l  - (Al-1))/(A-1)  ♦ CA1 

= ( A^  - ) ♦ o A1  = ( A^  ♦ A* ) 

which  contradicts  (3.16),  and  so  finishes  the  proof.  In  the  second  last 
equality,  we  used  the  fact  that  A-l  = (r2+l)/*2  - 1 = l/*2.  ||| 

Remark.  It  la  easy  to  see  that,  in  the  oorollary,  a^_^  ♦ a^,  ^ 0 
for  any  /u  in  case  the  nullspllne  s is  not  trivial.  Tor  if,  e.g.,  a_^  = 

*0  * 0.  than  s(r)  would  vanish  on  Ct-a^u(ltal),  s ^r'Vl1* 


i 


hence  e would  be  a polynomial  of  degree  < r on  that  Interval  and  vanish 
2(r-l)  times  there,  therefore  would  have  to  vanish  Identically  there. 

But  then,  we  would  have  s = 0 by  the  considerations  In  Section  2.  We 
can  therefore  conclude  from  the  corollary  that,  for  a nontrivial  null- 
spline  8, 

> const  A^ 

either  for  all  1 > 1 or  else  for  all  1 < -1,  with  a^  and  A as  In  the 
corollary  and  const  :=  j maxja^,  aQ\/(fc.  a)1  > 0. 

The  argument  for  this  corollary  would  have  been  simpler  bad  I been 
able  to  prove  that  every  p with  AP  = 0 can  be  written  as  a sum  p = p'+p" 
with  I1>0IPJ|2  < od  and  ^1<0lP^'l2  <®»  and  Ap'  = Ap"  =0. 

A minor  'variation  of  the  arguments  for  Lemma  3 and  Its  corollary 
allow  the  following  conclusion  of  Independent  Interest  In  the  study  of 
linear  difference  equations. 

Theor|m_^.  Let  A = (a^)  be  a bllnflnlte  matrix  which  represents 
a linear  map,  also  denoted  by  A,  on  ^(Z)  for  some  p e[l,  ®)  which  la 
bounded  and  bounded  below,  j.e..  there  exist  positive  E and  7 ao  that 

5,  Oallp  < IUa||p  < ^ flail p a C^(Z). 

2S  A is  a band  matrix,  l.e..  If 

■ :=  sup{|l-j|  s ajj  ^ o\  < ®, 

then  any  nontrivial  sequence  p for  which  Ap  = 0 must  Increase  exponent- 
ially either  for  Increasing  or  for  decreasing  1.  Explicitly,  there 
exist  an  Index  ft  and  a positive  eonst^  so  that,  either  for  all  1 > 
or  else  for  all  1 < /*., 

C P,p  * const. 

2ml< J<2m(l*l ) 3 

with 

A is  and  Kis  C/K. 

Thanks  are  due  to  Allan  Plnkus  for  questioning  the  necessity  of  an  addit- 
ional assumption  in  an  earlier  version  of  this  theorem. 
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J»j__SxEonentlal  decay  of  the  fundamental  spline.  Assume  that  the 

knot  sequence  Is  such  that  the  3.I.P.  is  correct,  l.e.,  has  exactly  one 

solution  s £nl  . for  every  atm(2 Z) . This  means  that  the  restriction 
a * »i 

map  R^,  when  restricted  to  Is  one-one,  onto,  and  clearly  hounded 

with  respect  to  the  sup-norm.  One  verifies  directly  (else  see  (4.2) 
below)  that  m$^  ^ Is  a closed  subspace  of  m(I),  hence  complete.  The  Open 
Mapping  Theorem  therefore  provides  the  conclusion  that  R^  Is  boundedly 
Invertible.  This  means  the  existence  of  some  const  so  that 

(4.1)  || SqIIjjj  < const  lal^  , all  ocm(Z). 

Let  k j.  be  the  1-th  B-spllne  of  order  k for  the  knot  sequ- 

ence t , normalized  so  that 

*!<*'  *■  w - t*i w"1'  - ^’i'1 

and  so  , comparing  with  the  ^splines  Introduced  at  the  beginning  of 
Section  3* 

S1 ,k ,i  = ((ti+k‘ti)/k)  Ml,k,i  • 

Prom  (3.3)*  or  already  from  [2], 

(4.2)  < <ViHA,  “ * *11  ptn(Z)* 

for  some  positive  constant  depending  only  on  k and  not  on  i . 

Since  (N^)^  la  a basis  for  . (In  the  sense  described  In  the  pre- 
ceding section).  It  follows  that  s ^ satisfies  s|^  * a If  and 
only  If  Its  3-spllne  coefficient  sequence  g satisfies 

(4.3)  EJHJ(tl)PJ  = °1'  a11  1 • 

while  s j Is  bounded  If  and  only  If  Its  corresponding  &-spllne  se- 
quence 3 Is  bounded,  by  (4.2).  We  oonclude  that  the  B.Z.P.  has  exactly 
one  solution  for  every  acm(Z)  Iff  the  matrix 

A :=  (lj(t&)) 

maps  faithfully  onto  4^.  We  oolleot  these  facts  In  the  following 
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Theorem  4.  The  bounded  Interpolation  problem  la  correct  If  and  only 


IT  the  matrix 


A = (N^)) 

provides  a faithful  linear  map  from  l (Z ) onto  ^ (Z ) . If  one  or  the  other 
of  these  conditions  holds,  then  A , being  trivially  bounded.  Is  bounded!? 
Invertible,  Since  A Is  also  a band  matrix,  of  band  width  m :=r-l,  It 
then  follows  from  Theorem  2 that  the  Inverse  of  A Is  also  given  by  a 
matrix,  (b^j)  say,  and  that 

|blj|  < const  all  1.3. 


with 

X:=(K/(l+tO)1/2m,  const  < KA2®,  «:  = QA"1^ 


since  = 1.  In  particular,  for  all  1 , the  function 


L1  != 


Is  then  a fundamental  spline  which  decays  exponentially  at  the  rate  X, 
and  the  solution  sQ  of  the  B.I.P.  for  given  a 6 m(2Z)  1b  given  by 


•a  = ZfjLi  • 


a aeries  which  converges  uniformly  on  compact  subsets  of  I. 


Ve  do  not  know  conditions  which  are  both  necessary  and  sufficient 
for  the  correctness  of  the  B.I.P.  . Since  correctness  implies  bounded- 
ness of  the  map  at— »sa#  we  obtain  from  [4;  Lemma  Of  Section  2]  the  nec- 
essary condition  that  the  local  mesh  ratio 

= sup  At. /At. 
i |i-  J ]=1  1 3 

be  finite.  If  the  local  mesh  ratio  is  Indeed  finite,  then  a simple  suf- 
ficient condition  for  uniqueness  of  the  interpolating  bounded  spline  is 
the  condition  that 


(4.4)  Is  •)  , 
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J, 


This  Is  connected  with  the  fact  that,  with  k = 2r,  the  r-th  derivative 
of  any  nontrivial  null  spline  grows  exponentially  In  at  least  one  direct- 
ion, as  described  In  Leona  3 and  Its  corollary.  Precisely,  we  have  the 
following 


Lemma  4.  If  m.  :=  sup  ht./ht.  <®,  and  there  exists  a bounded 

~ = |l-Jl=l  1 3 

nontrivial  null  spline  s In  ^ then  either  t^-r  -<d  or  ^ < ®. 

Proof.  Let  s = k be  the  nontrivial  bounded  nullspllne  in 

j.  Its  >th  derivative  Is  then  s^33  = with 


y(3) 

Ti 


( 


* 

(k-3)(r[3_1)  - 


l+k-  j V1 


i-t. ) , 


3=0 
3 >o 


This  Implies  the  estimate 


<4-5)  23  nax  W-31 

(see,  e.g.,  [4],  for  similar  considerations).  Write  now  the  r-th  deriv- 
ative of  s In  terms  of  the  somewhat  differently  normalised  B-spllnes 
:=  introduced  In  Seotlon  2, 


a 


(r) 


VA  * 


Then  ( ( t^-t^/r)1/2,  so  that,  from  (4.5), 

(4.6)  |p1l  < oonstr  lrfl-/(t1+r-t1)r"1/2  . 


By  the  corollary  to  Lemma  3 (in  Section  3),  we  may  assume,  without  loss 
of  generality,  the  existence  of  a positive  const  so  that,  with  m = r-1 , 

E |p, |2  > eonst  A3  , 3=2,3,... 

2m3<l<2m( 3+1)  x 

where  A:=  (1+I^)/k2  > 1 and  K<  D 2 , the  latter  a certain  constant  1&- 

• r 

dependent  of  i . In  con3unctlon  with  (4.6),  this  Implies  that 

oonat  A3  < constr  Y niax{(tl4r-t^)3,"2r  : 2m3<i<2m( 3*1)\ 

2 

£ eonstr>Y  (m^)r  ■ln{(t^+r-t^)3'“2r  s 2m3<l<2m(  J+1)\  , 
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where  we  have  used  the  fact  that 


"j-11*51  i i "I11'1'  • 

It  now  follows  that 

t1+r  - t ^ < const  A"1^21^  , l=2r,  2r+l,... 

and  therefore 

V * l2r  * ^I'lWIr  - *lr>  ' • * 111 

We  note  In  passing  that  the  argument  also  establishes  uniqueness 
In  case  either  t or  t^  Is  finite  as  long  as  the  local  mesh  ratio  Is 
for  some  <3  which  Is  greater  than  1 and  depends  on  A. 

We  are  now  ready  to  prove  Theorem  1. 

Proof  of  Theorem  1.  Slnoe  the  global  mesh  ratio  = 
sup,  ,ht,/ht,  Is  finite,  then.  In  particular,  I = (-«,®)  and  lemma  4 
Implies  that  R^  maps  ml^  ^ one-one  to  m(Z). 

Next,  we  prove  that,  for  each  1,  the  fundamental  spline  function 
Introduced  In  Lemma  1 decays  exponentially  away  from  t^,  l.e.,  for 
all  3 and  all  x«[tj,tj+^], 

(4.7)  |X.1(x)  | < const  )J1"3' 

for  some  const  depending  only  on  k and  M^,  and  some  Xt[0,l)  which  de- 
pends only  on  k.  It  suffices  to  consider  J > 1.  We  ha*e  L^(tQ)  = ® ^or 
n jt  1 , therefore 

Ll(x)  = (x”t^+l)*‘^x"tj*r^x,tJ+l"*,,t3+r^Ll 

= frtx.tj.! t'r>(t)dt/r!  . 

BL— 1 

By  Holder's  Inequality, 

« 

* (r/(t,„-x))1/2  Hlir>ll2,[Xlt^r]  > 

making  use  of  (3.3),  so  that,  from  the  oorollary  to  Lemma  2, 


-i\ - 


|L^( x)  | < constr  (EVh1^2)  con3tr  ||L^r^ !l 2 X'*"1 

with  X«  [0,1)  depending  only  on  k,  and 

E j=  supQ  &tn,  h :=  lnfn  HtR  . 

But  now,  from  Lemma  1, 

NLir^2  - constr  ^2/br  , 
and  (4.7)  follows. 

The  exponential  decay  of  all  fundamental  splines  L^  at  a rate 
which  does  not  depend  on  1 now  allows  us  to  construct  an  lnterpolant 
sa  In  ^ for  arbitrary  a €m (Z) , In  the  form 

8 a = £1°1L1  » 

which  satisfies 

II  sal^D  - oon3t  Hall® 

and  therefore  Is  In  m|^  ^ . | | | 

It  Is  clear  that  the  argument  for  Theorem  1 shows  the  existence  of 
a number  <?  > 1 (which  depends  on  k and  on  the  X of  Lemma  2)  so  that  the 
conclusions  of  Theorem  1 hold  even  If  we  only  know  that  the  local  mesh 
ratio  Is  less  than  A quick  analysis  of  the  constants  Involved  shows 
this  provable  ^ to  converge  to  1 very  fast  as  k Increases. 
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